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Abstract 

We resum distributions that are singular at partonic threshold (the 
elastic limit) in heavy quark production, in terms of logarithmic be- 
havior in moment space. The method may be applied to a variety of 
cross sections sensitive to the edge of phase space, including transverse 
momentum distributions. Beyond leading logarithm, dependence on 
the moment variable is controlled by a matrix renormalization group 
equation, reflecting the evolution of composite operators that repre- 
sent the color structure of the underlying hard scattering. At next- 
to-leading logarithmic accuracy, these evolution equations may be di- 
agonalized, and moment dependence in the cross section is a sum 
of exponentials. Beyond next-to-leading logarithm, resummation in- 
volves matrix-ordering. We give a detailed analysis for the case of 
heavy quark production by light quark annihilation and gluon fusion. 
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1 Introduction 

The factorization of cross sections into long- and short- distance functions 
gives QCD predictive power for a wide variety of inclusive, hard-scattering 
cross sections in hadronic collisions. There has been great progress in recent 
years in the calculation of perturbative corrections at next-to-leading order 
(NLO) for such processes, including Drell-Yan [|1|, direct photon 0, heavy 
quark and jet production 0]. Selected exact results at two or even more 
loops are also available, especially for e"'"e~ annihilation ^ and deeply inelas- 
tic scattering P|. The exploration of quantum chromodynamics, electroweak 
theory, and searches for extensions of the standard model, however, require 
information on yet higher orders. We need to develop methods to decide 
when higher-order perturbative terms are negligible and, when they are not, 
to estimate their importance. This motivates the systematic study of per- 
tubation theory at arbitrary order. One of the elements in this program is 
the resummation of singular distributions in factorized inclusive cross sec- 
tions. This resummation is most conveniently carried out in moment space. 
Closely related resummations organize large logarithmic corrections near zero 
transverse momentum for heavy systems such as Drell-Yan or heavy quark 
pairs. 

In this paper we shall give details of this resummation in hard-scattering 
cross sections induced by QCD itself, particularly heavy quark production, 
extending results for electroweak-induced hard scattering cross sections, such 
as the Drell-Yan production of vector bosons 0, ^, ||. Our starting point 
will be the factorization properties of these cross sections near threshold, 
and we shall make strong use of the close relation between factorization and 
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The resummations that we consider apply to processes of the form 

hiPi) + h^{p2) ^ T(g, y,x)+X, (1.1) 

where T{q,y,x) represents a system in the hadronic final state, produced 
with large invariant mass = Q and (collective) rapidity y, 



y=^lnr-^ . (1.2) 




Examples include the production of heavy quark pairs and high-p^ jets. The 
additional variable x represents the internal structure of the final state, for 
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example the angle 6 between the quark, or jet, direction in the final state 
and the beam axis. Such cross sections may be written in factorized form as 



h2 



dQ"^ dx dy 



J dXadXb 0a//ii (Xa, /i^)0b/^2 i^b, yW^) 




Q 



c.,2/-2/p,X, — , (1.3) 

where S is the center-of-mass (cm.) energy squared of the hadrons, and 
where the 0's are parton distributions defined in an appropriate scheme, and 
evaluated at factorization scale /i. The function Hab is the hard-scattering 
partonic cross section (the "hard part") that describes the production of 
the final-state system T in the collision of partons a and 6, which carry 
momentum fractions Xa and x^ of hadrons hi and /;,2, respectively. Rapidity- 
dependence in Hab is through the combination y — i/p, where yp is the partonic 
rapidity in the hadronic cm. frame, 

yp = l^-(-)- (1-4) 



2 \xb 

The hard parts Hab in Eq- ( |1.3| ) are calculable in perturbation theory. In 
general, they are not smooth functions of their arguments, but combinations 
of functions and singular distributions. The reason for this somewhat compli- 
cated structure is that they summarize cancellations between cross sections 
with gluon emission and with virtual gluon corrections. This cancellation 
(in QCD as in QED) is only manifest at the level of sufficiently inclusive 
cross sections. For factorized cross sections such as Eq. (p.. 3D, this require- 
ment manifests itself in the occurrence of ("plus") distributions, which are 
singular for z = 1, where 

z = ^, (1.5) 

s 

with s = XaXbS the invariant mass squared of the partons that initiate the 
hard scattering. We shall refer to z = 1 as "partonic threshold" Q or some- 
times as the "elastic limit". 



-'^We emphasize that by partonic threshold, we refer to the cm. total energy of the 
incoming partons for a fixed final state; heavy quarks, for example, are not necessarily 
produced at rest. 
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Singular distributions are most easily organized in terms of moments with 
respect to r = Q'^/S. As we review below, for the Drell-Yan and related 
cross sections logarithmic dependence on the moment variable exponentiates 
This moment must be inverted to derive the physical cross sec- 
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tion. There are a number of proposals on the best way to carry out this 
transformation, or otherwise use the resummation, which have been put for- 
ward in the context of Drell-Yan [12, 13, O] and heavy quark production 



|I5| , |1^, |1^, |1^, |1^, Our analysis here will not deal directly with this 
question; our results will remain in moment space. In this paper, we shall 
be concerned with the extension of resummation to cross sections that are 
initiated by QCD itself. Our explicit applications will be to resummation 
at partonic threshold for heavy quark production, but the method is much 
more general. A straightforward application of the method, which involves 
only the replacement of the moments by Fourier transformation to impact 
parameter space, is to the overall transverse momentum of the heavy quark 
pair, again extending results for Drell-Yan pairs and hadrons in e"'"e~ 
annihilation [p^ . 

It has been recognized for some time that leading logarithms in moment 
(or impact parameter) space exponentiate for heavy quark production |]T3 



1^ in the same manner as for Drell-Yan, reflecting the factorizability 

Beyond the 



of leading singular distributions from the hard subprocess 
level of leading logarithms, however, singular distributions in QCD-induced 
cross sections depend coherently on color exchange in the underlying hard 
scattering. The interplay between the color structures of short and long 
distance interactions will be a unifying thread of the discussion that follows. 

The next section starts with a brief review of the factorization of the in- 
clusive Drell Yan cross sections near partonic threshold. We shall describe 
the resummation of the singular distributions in these cross sections starting 
from their factorization properties. We also observe the field-theoretic con- 
tent of the functions into which the cross section factorizes, relating them to 

26|, 13] 



23, 24, 25 



We then ex- 



matrix elements of fields and of Wilson lines 
tend the factorization, and hence resummation, formalism to QCD-induced 
inclusive cross sections. The possibility for, and the broad outlines of, such 
an extension has recently been discussed [10|; here, the practical formalism, 
and its field theoretic basis, are discussed in some detail, expanding on a 
summary given in Ref. ||2^ (see also [^]). To organize sensitivity to color 
exchange at short distances beyond leading order in ln(l — z) (or in loga- 
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rithms of the corresponding moment variable), we derive a matrix evolution 
equation that controls all logarithmic enhancements at threshold. At next- 
to-leading logarithm (NLL) in the exponent, the matrix equation may be 
diagonalized in a manner reminiscent of leading order evolution in deeply 
inelastic scattering. At NLL, the cross section is a sum of exponentials, each 
representing an eigenvector of the anomalous dimension matrix. Beyond this 
level of accuracy, however, the resummation can only be carried out in terms 
of ordered exponentials. Resummation at NLL requires the computation of 
the anomalous dimension matrix for each partonic subprocess, such as light 
quark annihilation into heavy quarks or gluon fusion into heavy quarks. 

The third section describes the calculation of the anomalous dimension 
matrix for light quark annihilation into heavy quarks, supplementing the 



brief description of Ref. p8[. Section 4 extends these results to the very 
important case of gluon fusion. The paper concludes with a summary and 
the outlook for further progress. 



2 Factorization and Resummation 

2.1 Resummation for Color Singlet Hard Scattering 

The resummation of singular distributions for the Drell-Yan, and related 
electroweak-induced hard scattering processes has been understood for some 
time 0, H, If. In this subsection, we review the basic results, and elaborate 
somewhat on the operators in QCD that play a role in the underlying factor- 
ization. For our purposes, it will be useful to generalize the standard results 
slightly, to the production of an arbitrary color-singlet system, denoted B 
below, of invariant mass squared Q^. We assume in particular that B may 
be produced by gluon fusion as well as quark annihilation. An important 
example is Higgs production [0. For the purposes of the following discus- 
sion, however, we shall refer to this class of production processes as Drell-Yan 
cross sections. 

Resummation is most readily formulated in terms of moments of the 
inclusive Drell-Yan cross section da/dQ"^, with respect to the variable r = 



6 



xuu^B {n, Q/fi, asifi^)) + Oil/N) . 

(2.1) 

Here, ej is the fractional quark charge, do (= 47ra^/(9(5^5') for the standard 
Drell-Yan cross section) is a Born cross section, and the moment functions 
are ^ 

4){N,i2^)= f dxx^^^(j){x,fi^), (2.2) 
Jo 

and similarly for the hard-scattering function U{i_^^{z, Q/ fi, ag^fi^)) in terms 
of the variable z. In Eq. ( |2.1| ) we have neglected hard-scattering functions 
associated with the gluon-quark combination, with two quarks or antiquarks, 
and with quarks and antiquarks of different flavors, since these do not give 
rise to singular distributions at 2; = 1. All such corrections, which are smooth 
at z = 1, contribute to the moments only by corrections that vanish as powers 
of for — > 00, as indicated. Singular distributions, however, give rise to 
logarithmic A^-dependence, as in 

. AT + . . . (2.3) 



f'dz 


"ln'"(l-z)" 


1 - z 


Jo 



m + 1 

Thus, we should think of the large- A^ behavior as a diagnostic for singular 
distributions in 1 — Physical cross sections, of course, always involve the 
inverse transform of the moments. 

The functions tUj^^^ (2;, QZ/i, as(/i^)) whose moments appear in Eq. (|2.1| ) 
are normalized to 

uJu^B = 5(1 -z) + 0{as) . (2.4) 

We shall see that the logarithmic A^-dependence of moments of u exponen- 
tiates into a useful form. 

Before sketching the physical basis and derivation of this result, we ob- 
serve that, as shown in the same functions ojfi^s as in Eq. ( p.l|) contain 



the singular distributions at z = 1 for fixed rapidity 
dahihi^B J f dxb 




dQ'^dy .^-^ Jt J Xa Xb 

2 Xb 



(2.5) 
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Thus, we may compute singular threshold corrections for the cross section at 
fixed rapidity from the inclusive cross section. 

The resummation of singular distributions in 1 — z may be carried out in 
a number of ways, but all are based on the observation that the perturba- 
tive function ojfi may be computed in dimensionally regulated perturbation 
theory by choosing the external hadrons to be partons. For us the most infor- 
mative approach is based on the ability to reexpress moments of the inclusive 
Drell-Yan cross section Eq. (|2.1|), with, for instance, /ii = g, /12 = ^ [05 i]- 
In this form, all logarithmic behavior is factorized into alternate parton 
distributions if), along with a "soft-parton" function f/, which describes the 
exchange and emission of soft gluons by the annihilating quark pair 0, 

+ 0{l/N) . (2.6) 

This factorization is illustrated in Fig. 1. i7 is an infrared-safe function that 
is independent of N. Let us discuss the remaining functions, which absorb 
the A^-dependence, in turn. 

The modified parton distributions ip in Eq. ( p.6| ) are conveniently defined 
in the partonic center of mass system at fixed energy, rather than a light-like 
momentum fraction. For i = q, for example, we have0 

1 r°° ■ ^1 

^P,/,ix, 2po//i, e) = J dy^ e-™°(g(p)|g(yo, 0) ■ 7 g(0)|g(p)) , 

(2.7) 

where the matrix element is evaluated in = gauge in the partonic cm. 
frame. The vector n is lightlike in the opposite direction from p^, so that 
for p in the ±3 direction, n ■ 7 = 7'''. ipg/g and ipg/g may be constructed 
similarly. Defined in this fashion, ip absorbs all the collinear divergences of 
the Drell-Yan cross section in the r — >■ 1 limit 0. Its moments may thus be 
factorized into a product of moments of the corresponding light-cone parton 
distribution 0, defined in any scheme, times an infrared safe function. 

^This form includes an explicit spin average. Single-particle states are normalized by 
(0|g(0)|p, s) = u{p, s). The factor of 1/2^^^ then normalizes the lowest order to 6(1 — x). 
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Fig. 1. (a) Factorization for the Drell-Yan process near partonic threshold, 
Eq. ( |2.6|) . (b) The soft-gluon function U. The double lines represent ordered 
exponentials in Eq. (|2.10|) . 
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The function [7*^*^ represents the couphng of soft gluons to the annihilating 
partons. So long as these soft gluons are not collinear to the active partons, 
this couphng is well approximated by ordered exponentials, or Wilson lines, 
which we denote as 

$^(A2, Ai; x) = P exp (^-ig £^ dX'(3 ■ A{X'(3 + x)^ , (2.8) 

with A the gauge field in the appropriate representation of the gauge group, 
and (3 a velocity four-vector, in the direction of the particle whose interac- 
tions with soft gluons are being approximated. In the case of a quark, the 
gauge fields are in the fundamental representation, for a gluon in the adjoint. 
The operator P orders group products in the same sense as ordering in the 
integration variable A, with the A's with lower values of A to the right. In 
this notation, $/3(A2, Ai; x) represents (for example) a quark that propagates 
from X + XiP to a; + X2P, while $/3(Ai, A2; x) represents an antiquark with 
the same propagation. The construction of U^^^ from the $'s is actually a 
two-step process, which we now describe. 

Starting with the Wilson lines, we first construct a (dimensionless) "eikonal 
singlet cross section" 

a^^? ((1 - u;)Q//i, asifx'), e) = ^ | dxo e'(i-)«-« (0|f «,,(xo, 0)" 

xT[W,,,M]\0), (2.9) 

where W{x) is the product of Wilson hues representing the pair, coupled by 
a color-singlet vertex at the point x, 

WbiMi^) = ^aia2 '^'/32(-oo,0;a;)b2a2$/3i(0, -cx);x)aif,, • (2-10) 

In Eq. ( ^.91 ), T represents time-ordering, and T anti-time ordering. When 
dealing with QCD hard-scattering below, we shall generalize this product, 
by incorporating nontrivial color structure at the vertex between the lines. 
The quantity {1 — w)Q is the energy emitted by the Wilson lines [27, 32 



which is the relevant measure of phase space near partonic threshold in the 
center of mass frame of the colliding partons |^. In Eq. ( p.9| ), an appropriate 
normalization is understood, so that at lowest order ctj^y'* = S{1 — w). 

The eikonal cross section represents the emission of soft gluons from the 
annihilating pair. The scale invariance of the Wilson lines implies, however. 
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that it contains ultraviolet, in addition to infrared and collinear, divergences, 
corresponding to the renormalization of the composite operator introduced 
by the vertex in Eq. ( |2.1(J| ) where the eikonal lines join [^, |2^, |2^. A 
natural "scheme" for this renormalization is specified by demanding that the 



integral of (J^y\^ ~ from w = to w = 1 he normalized to unity. 



dw cr 



(cik) 
DY 



1 - w)Q/jJ,, OsifM^ 



1. 



(2.11) 



That is, we demand that all corrections to cr^Y ^ vanish for the fully inclusive 
cross section. With this normalization, beyond lowest order cr^y^ is a sum 
of "plus" distributions in 1 — w. 



5{l-w) 

fi,as{fJ^ ),e 



n>0 



ln"(l -w) 
1 — w 



(2.12) 



where the coefficients s„ are power series in the coupling. 

Even after renormalization, cty^y^ includes the couplings of fast-moving 
gluons that are collinear to the incoming Wilson lines. Such gluons, whether 
virtual or real, are already included in the distributions ip, and their con- 
tributions should be removed from the soft gluon function, f/^*\ to avoid 
double counting. This is possible because noncoUinear soft gluons factor 
from collinear gluons through the use of Ward identities ||3^. Physically, soft 
radiation cannot distinguish between physical and eikonal light-like sources. 
Thus, moments of factorizable in the same manner as the quark- 



antiquark cross section cross section, Eq. (2^), into moments of the (identi- 
cal) soft-parton function U times moments of "jet" functions ji_2, analogous 
to the -^'s, Eq. (|2.7|), but with the incoming partons represented by products 
of eikonal lines (in Aq = gauge) as 



Ji 



Wi)Q 



^ r dyo e-^^i-'"^)^^'' 

x(0|Tr|f[4^(0,-oo;y)] 
xr[$^,(0,-oo;0)] 1 |0), 



(2.13) 
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with = {no, 0) a vector at the spatial origin. The "antiquark" distribution 
j2 is defined analogously, but with the Wilson line in the P2 direction and 
ordered from to —00. 

In these terms, the eikonal version of the moment relation Eq. ( |2.6| ) is 



;^{eik) 



a 



DY 



[7» 



'Q_ 



Q_ 



32 



'Q_ 

Nix 



a^(/i^),e 



(2.14) 



which may be used to construct the soft function U . This prescription elim- 
inates diagrams that have divergences from gluons collinear to the external 
eikonal lines of ?7, a fact which we shall use below. Because the jVs absorb 
all collinear divergences, V is an infrared safe sum of plus distributions. 
Logarithms of N exponentiate in each of the factors of Eq. (|2.6|) . As 



emphasized in Ref . ||10[ , this is a result of the constraint that the product is 
independent of both the gauge choice and the factorization scale, although 
its individual factors are not. Such combined constraints organize logarithms 
associated with both soft and collinear momentum configurations, and are 
at the basis of Sudakov resummation 0]. Comparing Eq. (|2.1|) with h\ = i, 
h2 = i with Eq. ( |2.6| ), we can solve for Co{N), in terms of the calculable 
A^-dependence oi ip, (J) and U^^^ 0, 



a;, 



{N) 



i'i/iiN,l,e) 
4>^/^{N,Q^e) 



(2.15) 



where we have used ipi/i = ipi/j, and similarly for the 0's, and have suppressed 
various arguments. 

The result of this procedure applied to the Drell-Yan cross section is 



^ii-.BiN) = A{asiQ^)) exp 



— I dz 
Jo 



1 ^N-l 



(1-^)2 A 



(2.16) 
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where for DIS (MS) factorization scheme ms = 1 (0). A{as) is a constant 
function of the couphng, and the gi^2 are finite functions of their arguments. 
g["'''^ and (72"'''' are universal among hard cross sections and color structures 
for given incoming partons a and b, but depend on whether these partons 
are quarks or gluons. To reach the accuracy of next-to-leading logarithms in 
the exponents, we need gi only to two loops, with leading logarithms coming 
entirely from its one-loop approximation, and g2 only to a single loop. More 
explicitly, we take |jl3[ 



(2.17) 



g^r'^ = (a + a 




with Ci = Cf (Ca) for an incoming quark (gluon), and with K given by 



K 




9 ^' 



(2.18) 



where n/ is the number of quark flavors. g2 is given for quarks in the DIS 
scheme by 

(2.19) 



and it vanishes in the MS scheme. The resummed A^-dependence thus de- 
pends on the factorization scheme, and is quite different between DIS and MS 
schemes p, |Ty|. This difference is in principle compensated for by differences 
in the A^-dependence of the distributions themselves. 

The leading logarithms in the exponent of Eq. ( |2.16[ ) are from the ratio 
of 'ip to (p in Eq. ( p.lSj) . The function f/^*) contributes at next-to-leading 
logarithm. To achieve this accuracy we need only a one-loop calculation, 
and we find M 



exp 



exp 



dz — 
1 

1 z 

dz — 
1 



N-1 



z/«(a.((l-z)^g^)) 



z 



Z 



t''^^^' ^/3(^?(Ag2))z.»[a.(Ag2)] 

'(1-2)2 X 



W^)\a,{{\-zr^Q')\ 



(2.20) 



with v^'^'^^ = 2Cf,a{ois/'^) + • • • In the second form, (3{g) is the QCD beta 
function, and we have followed H in reexpressing ^Z*^*-* into the form of Eq. 
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( |2.16| ). In this connection, we note that in axial gauges receives noncan- 
celhng contributions only from the "gauge" term k^^ /{n ■ kY in the gluon 
propagator. At least to this order, it is natural to use Ward identities to 
factor such terms into the parton distributions if). 



2.2 Resummation for QCD Hard Scattering 



Generalizing Eq. ( ^.51 ), we consider the production of a pair of heavy quarks 
(or jets) at rapidity y, and at scattering angle 9 in the pair center of mass 
frame, as 



dQ'^dyd COS 9 



f=q,q,g •' ^ ^" 




5[z -; ) 5 [y — 



Xq^xi)S I \ z Xb 



+ .... (2.21) 

Corrections, due to other partonic flavor combinations, are less singular by 
a power of 1 — 2 as 2 1 in the hard scattering, and their moments are 
suppressed by powers of 1/A^, as in Eq. (|2.1| ). For definiteness, we shall 
concentrate on heavy quark production, but our considerations are more 
general. As in the case of Drell-Yan cross sections, the singular behavior of 
the hard-scattering functions a (the analog of Ugg above) at fixed rapidity may 
be computed from the inclusive partonic cross section Following Ref. 
[0, we begin with the factorization properties of heavy quark production 
and other QCD hard-scattering processes in the elastic limit. We shall see 
that these depend upon the directions of the outgoing heavy quarks in the 
partonic center of mass frame. As a result, resummation must be carried out 
at fixed directions for the heavy quarks. The reasoning is then very similar 
to that for the Drell-Yan process discussed above. Note in particular, that 
the partonic subreaction remains an annihilation (or gluon fusion) process. 

To apply our methods to heavy quark and other QCD production mech- 
anisms, Eq. ( p.6| ) must be generalized in a number of ways. First, the un- 
derlying hard process involves not only a singlet color configuration, but any 
color tensor that may be constructed from the color representations of the 
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incoming partons. Second, the outgoing heavy quarks (or other particles), 
as well as the incoming light partons, act as sources of soft gluons. Our 
analysis is simplified, however, by the following observation. Near partonic 
threshold, the emission of light quarks into the final state is suppressed by 
factors of 1 — aside from pairs produced by soft gluons. As a result, al- 
though different color structures may "mix" due to soft gluon emission, the 
flavor structure of the hard scattering is unchanged by all corrections that 
are singular in 1 — 2;. Thus in the case of heavy quark production, we may 
analyze light quark annihilation and gluon fusion separately. 

To take these new features into account, we generalize the Drell-Yan 
factorization to QCD hard scattering initiated by light partons / and / with 
the following expressions, which are the analogs of Eqs. (for hi = f, 

h = f) and dp). 



'rfrr--^^^(r,Q^e) 
clQ^^ 

= (j)f/f{N, /i^ e) %-,qq(A^, <5//i, Q, as(/i^)) <Pf/f{N, /i^ e) 

= J2K ( ?)«s(/i^) ) S,ji(-^,as{fi'^)] hi ( ^,a^(/i^^ ' 
ij " 



fi J \ fiN ) \fi 



x^f/f i^N, ^, ej l^N, ^, ej + 0{1/N) . (2.22) 

This factorization is illustrated in Fig. 2, for the process 

fiPa) + fiPb) -> QiPl) + Q{P2) , (2.23) 

in which // represents a pair of light quarks or gluons that annihilate or 
fuse, respectively, into a heavy quark pair. 

This factorization differs from Eq. ( |2.6| ) only in taking into account color 
exchange at the hard scatterings, hi in the amplitude, and h*j in the complex 
conjugate. For simplicity, we omit the partonic indices /, / in hi and /i}, 
although these functions certainly depend on incoming flavor. For quark- 
antiquark annihilation, the indices I and J range over two values, corre- 
sponding to annihilation in a singlet or octet in the s-channel (or t channel). 
For other hard-scatterings, such as gluon-gluon, the same pattern applies, 
but with I and J ranging over color tensors appropriate to possible com- 
binations of gluonic representations in the initial state with (heavy) quark 
representations in the final state. 



15 




Fig. 2. (a) Factorization for heavy quark production near partonic threshold, 
Eq. ( p.22| ). (b) The soft-gluon function Sjj, in which the vertices T hnk 
ordered exponentials as in Eq. ( p. 24]) . 
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The color tensor of the hard scattering serves to tie together the ordered 
exponentials to which soft gluons couple, generalizing the singlet vertex join- 
ing the eikonal lines in Eq. ( p^.lOD above, 



x{Ti 



^ J di,da,djdi 



^Va (-00, 0; X)c,^da^^^{0, -00; X)d,, 



(2.24) 



where is in the color representation flavor /, with velocity Va- Notice that 
only the Wilson lines of the incoming partons have light-like velocity vectors. 
Next, by analogy to Eq. ( [^.9D for the eikonal cross section in Drell-Yan, we 
define. 



a 



(eik) 
JI 



[1 - w)Q/fi, asijJ,' 



I Jdxoe 
xT 



(oif [^^(xo,o)|,^; 

|0) . (2.25) 



As in the Drell-Yan case, this cross section includes unphysical collinear 
divergences associated with its incoming, lightlike eikonal lines. These are 
eliminated in the soft function Sjj by refactorizing aff^^ in the same manner 

as (T 7 

Eq. (|2.14| ), through moments. 



~(eik) 
JI 



a 



{N)=Sjj{N)~niN) UN). 



(2.26) 



Now, comparing Eqs. ( |2.6| ) and ( p.22| ), we see that the A^-dependence 
of the cross section for heavy-quark production from a set of initial-state 
partons is given by the following generalization of Eq. ( p^.lS] ) , 



(AT) 



^f/fiN,l,e) 



1 2 



h*j{Q/fi) Sjj{Q/Nfi) hj{Q/fi) , (2.27) 

IJ 



where the sum over / and J is, as usual, over color tensors. To determine the 
singular behavior at partonic threshold, we only need to resum logarithmic 
A^-dependence in each Sji. 

The soft matrix Sji depends on A^ through the ratio Q/{N^). Its A^- 
dependence, then, can be resummed by renormalization group analysis ^2 



|2^ , |25| , p^ , 0. The composite operator in Eq. ( p.24| ), which defines Sji, 
requires renormalization. Nevertheless, the product h*jSjihi of the soft func- 
tion and the hard factors needs no overall renormalization, because the UV 
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divergences of Sji associated with the new composite operator are balanced 



by those of the /I's. This situation is summarized by the conditions |T0 

1=1 

i=l 

^fi = {Zs)jbSbaZs,ai, (2.28) 

where h^'^^ and S^^^ denote the unrenormahzed quantities, is the renormal- 
ization constant of the ith incoming partonic field external to hj and Zsjj 
is a matrix of renormalization constants, which describe the renormalization 
of the soft function, including mixing of color structures. Zs is defined to 
include the wave function renormalization necessary for the outgoing eikonal 
lines that represent the heavy quarks. 

From Eq. ( p.28| ), the soft function Sjj satisfies the renormalization group 
equation ||2^ 



+ ^^^^^) ^ '(^s)jbSbi - Sja{Ts)ai , (2.29) 

where the anomalous dimension matrices may be found by explicit renormal- 
ization of the soft function. Following [ES, we find it convenient to compute 



the matrix in a minimal subtraction scheme. To be explicit, taking e = 4 — n 
the matrix of anomalous dimensions may be found at one loop from the 
matrix of renormalization constants by 

Tsig) = -|^Res,^o^s(^, e). (2.30) 

The determination of the anomalous dimension matrices will be the sub- 
ject of the following sections, for heavy quark production through the anni- 
hilation of light quarks and gluon fusion. 

The solution to the renormalization group equation (|2.29|) organizes log- 
arithms of in the soft function. In general, this matrix equation involves a 
scale-dependent mixing of tensors, and has no closed expression. At the level 
of leading logarithms of N in Sji, and therefore at next-to- leading logarithm 
of N in the cross section as a whole, we may simplify by choosing a color 
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basis in which the anomalous dimension matrix is diagonal, with eigenvalues 
A/ for each basis color tensor labelled by /. At this level of approximation, 
we have 



S 



ji 



S 



JI 



Q 

\Nfii 
X exp 



Ml 1^ 

(2.31) 



The coherent effects of soft gluon emission are all contained in this expression, 
which we now employ to determine ctjj^qq, Eq. ( p.27| ), to NLL accuracy. 

Our first step is to rewrite Eq. ( ^.27| ) slightly, to isolate differences and 
similarities to the Drell- Yan result of Eq. (|2.15|) , 



ff- 



<Pf/f{N,Q^,e)\ 
^ XuKiQ/f^) Sji{Q/Nfi) hr{Q/fi) 
U(n(Q/Nfx) 



(2.32) 



We have simply multiplied and divided by U^^^{Q/Nfi), defined in Eq. (|2.14|) . 
The first factor in Eq. ( |2.32| ) is just the resummed singlet cross section ( p.l6| ). 
The new factor is to be treated using Eq. (|2.31|) and the explicit expression 
off/, Eq. (^. 

Combining the results just discussed, and picking a color basis in which 
moments with respect to z exponentiate, we find, at NLL accuracy 



a 



= A'ias{Q'))h}{l,as{Q')) Sji(l,as{[Q/N]')) hi{l,asiQ' 



X exp 



Ejj{N, 9,Q' 



(2.33) 



where A' is an overall constant and where the color-dependent exponents are 
given by 



Ejj{N,9,Q' 



dz— 



N-l 



- 1 



A 
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(2.34) 

g^^ summarizes soft logarithms that depend directly on color exchange in the 
hard scattering, and hence also on the identities and relative directions of the 
colliding partons (through 6*), both incoming and outgoing. As pointed out 
above, one-loop contributions to g^ may always be absorbed into the one- 
loop contribution to g2 and the two- loop contribution to gi 0. Because 
g\l^ depends upon /, however, it is advantageous to keep this nonfactoring 
color-dependence separate. 

Given a choice of incoming and outgoing partons, next-to-leading loga- 
rithms in the moment variable exponentiate as in ( |2.33| ) in the color tensor 



basis that diagonalizes T^'g^J^^''^ , with eigenvalues A/. The resulting soft gluon 
anomalous dimension g^^ is then simply 

g\l\a,,e] = -h[as,e] + ]^v^'\a,) , (2.35) 

where the eigenvalues are complex in general, and depend on the relative 
directions of the incoming and outgoing partons, as shown. 



3 Applications to qq QQ 

These considerations may be illustrated by heavy quark production through 
light quark annihilation, 

qiPa) + qiPb) QiPi) + QiP2) . (3.1) 



In this case, as in quark-quark elastic scattering p4| , Pq] , the anomalous 
dimension matrix is two-dimensional. Following Ref. |3^ we define the in- 
variants 

s = {Pa+Pbf, h = (pa - Pif - m'^ , ui = {pb - pif - rri^ , (3.2) 
with m the heavy quark mass, which satisfy 



s + ti + ui = Q 



(3.3) 
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at partonic threshold. We also define dimensionless vectors ff by 

= Q^r (3.4) 

which obey vf = for the light incoming quarks and vf = / for the 
outgoing heavy quarks. Note that Q satisfies the kinematic relation s = 2Q^. 

We now calculate the matrix Ts{g) for this case. Renormalization asso- 
ciated with incoming light partons is factorized into the wave functions ip of 
the previous section. We therefore need consider only diagrams which are 
not part of the eikonal "jets" ji, as in Eq. ( p.l3|) . In particular, we may omit 



the UV structure of the incoming light-like eikonal lines. 

The UV divergent 0{as) contribution to each Sji is the sum of graphs 
illustrated by Fig. 3. Most interesting are the "vertex corrections" in Fig. 
3(a). The counterterms for S are the ultraviolet divergent coefficients times 
our basis color tensors, denoted q 

S2 = CiZs,12 + C2^5,22 , (3.5) 

where we will specify a convenient basis for the Cj below. 

In our calculations we employ the general axial gauge gluon propagator, 

D^\k) = —^N^'^{k), iV^-(fc)=^^-- ^ +n^i^,, (3.6) 

k'^ + le n ■ k [n ■ ky 

with the gauge vector, and eikonal rules for all external lines (Fig. 4). We 
now review the notation of Ref. which enables us to summarize results 
for the diagrams of Fig. 3 in a unified manner. 

If we denote a typical one-loop correction to c/ as uj^^\5iVi^ 5jVj, Aj, Aj), 
where 6 = +1 (—1) for the gluon momenta fiowing in the same (opposite) 
direction to the momentum of Vi and A = +1 (— 1) for Vi corresponding to a 
quark (antiquark), then we have: 

J (27r)" q'^ + ie \ {6iVi ■ q + ie){5jVj ■ q + ie) 



AiVi ■ n P A,f ■ n P 9 P 

+ n 



{6iVi ■ q + ie){n- q) {5jVj ■ q + ie) {n ■ q) {n ■ q) 



2 



(3.7) 
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(b) 

Fig. 3. One-loop corrections to Sji for heavy quark production through 
hght quark annihilation: (a) vertex corrections; (b) heavy quark self-energy 
graphs. 
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v7(-v-k+i8) -v7(v-k+i8) 




-v7(-v • k+i 8) v7(v • k+i 8) 

Fig. 4. Feynman rules for eikonal lines in qq — > QQ. The gluon momentum 
flows out of the eikonal lines. Group matrices at the vertices are the same as 
for quark lines. 



where c/ is a color tensor. P stands for principal value, 



1 



1 



(g ■ nY 2 \{q ■ n + leY 



1 



(— g ■ n + leY 



We rewrite ( 



as 



Ii{6iVi,6jVj) - ^l2{SiVi,n) - ^hiSiVi, -n) 
']^h{5jVj, n) - ]^h{5jVj, -n) + U{n^) 



where the overall sign is given by 

Sij = Aj Aj 5i 5j. 



(3.8) 



(3.9) 



(3.10) 



We now evaluate the ultraviolet poles of the integrals. For the integrals 
when both Vi and Vj refer to massive quarks we have (with e = 4 — ra) 



jUV pole 






TT e 


jUV pole 
-'2 








TT 


e 
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rUV pole 



jUV pole 
-'4 



_ T 

TT e 



(3.11) 



where the Lp is the famihar velocity-dependent eikonal function 



2mVs / 1 
- m 



l + f3 



+ ni 



(3.12) 



with /3 = Y^l — A.W? I s. The Lj and are rather complicated functions of 
the gauge vector n. We will see shortly that their contributions are cancelled 
by the inclusion of self energies. Their explicit expressions are: 



-n) 



Li{-n)] , 



(3.13) 



where 

Li{±n) 



Vi ■ n\ 



2 ^ (vi ■ nY — 2rn?'n? /s 

5(±n) 2rn? / s — \vi ■ n\ 

In ' 



[vi ■ nY — Irn^rfi j s 



+ In 



S{±n) 2w? /s — \vi ■ n\ + ^J {vi ■ nf — 2rf?ri? j s 

(5(±n) — \vi ■ n\ 
5(±n) n2 - 



(w; • n)^ — 2m'^n'^ / s 
n\ + ^Jiyi • nY — 2rTn?'n?/Sj 



(3.14) 



with 5{n) = \vi ■ n\/{vi ■ n). 

When Vi refers to a massive quark and Vj to a massless quark we have 



UV pole 



I, 



jUV pole 



jUV pole 
-'3 



jUV pole 
-'4 



2^ 



7 + In 



2m? 



ln(47r) 



TT e 



^||_l[^ + ln(^^.)_ln(47r)]}, 



1 



(3.15) 
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where 



(3.16) 



and Vij = Vi ■ Vj. Note that the double poles cancel in the sum over the J's 
in Eq. (pj). 

Finally, when both Vi and vj refer to massless quarks we have 

1 



jUV pole 
-'1 — 


as ^ 

IT 


' 2 


tUV pole 
-'2 — 


as 
27r ^ 




jUV pole 

-'3 — 


as 
27r ^ 


r2 


jUV pole 
-'4 — 


a. 


1 




vr 


e 



7 + In 



ln(47r) 



--[7 + ln(z/,)-ln(47r)]| 



--[7 + ln(i.,)-ln(47r)]| 



(3.17) 



Again, note that the double poles cancel. 

Our calculations are most easily carried out in a color tensor basis con- 
sisting of singlet exchange in the s and u channels, 

Cl = SabSu , C2 = Sa2Sbl ■ (3.18) 

The color indices for the incoming (light) quark and antiquark are a and 
b, respectively, and for the outgoing (heavy) quark and antiquark 2 and 1, 
respectively. 

In the basis ( ^.181 ) we find from Eq. ( |2.30D that the contributions of the 
diagrams in Fig. 3(a) to the anomalous dimensions are 



r5,ii 
r5,2i 

rs,i2 



a 

TT 



Vab 



In — - L/3 - - \n{i'ai^b) - Li- L2 



5,21, 



as_ 
27r 
as 
2^ 

as_ 

TT 



In 



Va2Vbl\ 
ValVb2) 
Vab 



In ( — ) - \n{vaiVb2) - Lp + \n. 



2m2 



Til 



\n{va2Vbi) - In 



2m2 



\n{iJai^b) - Li - L2 + 2 



N 



S,12- 



(3.19) 
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The matrix depends, as expected, on the directions of the Wilson hnes, which 
may be reexpressed in terms of ratios of kinematic invariants for the partonic 
scattering. It is also gauge-dependent. 

We now combine this matrix with the (color-diagonal) contributions from 
the heavy-quark self-energy graphs in Fig. 3(b). The contribution of the self- 
energy graphs is 

(3.20) 



a 



— Cf{Li + L2 - 2) 61J 



which cancels the gauge- dependent Lj of Eq. ( p.l3| ). In addition, recalling 
Eq. ( |2.35| ), we subtract half of the anomalous dimension 
associated with the soft gluon function U^^^ of the Drell-Yan cross section 
(see Eq. ( p.32|) ). This shifts our matrix by an additional — {as 11^)0 pSij. 
Combining these contributions with Vs of Eq. (|3.19|) , the full soft-gluon 
anomalous dimension matrix becomes 



J- 5,11 

r5,2i 
rs,i2 
rs,22 



7r 



In 



Vab 

2 



■ ln(z/aZ/b) 



TXl 



-^r5,2i, 



2^ 
27r 



^In 
In 



Va2Vbl 
ValVb2 
Vab 

2 



- \ia{vaiVb2) - Lfs + In 



2m2 



— TCZ 



TT 



\n{va2Vbi) - In 



2m2 



-^r5,i2. (3.21) 

In the = gauge we have I'a = ^b - 
mass frame. Using the relations Vab = 1, V12 
and Vib = V2a = —Ui/s, the anomalous dimension matrix becomes in terms 
of our invariants s, ti, and Ui : 



1/2 in the partonic center of 
1 - 2m^/s, Vai = Vb2 = -ti/s 
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r5,i2 



In 



TT 



A. 



In 5 



+ 



Ca 



-In 



' rr?s' 

IT, 



Lp - ni 
-Lp-l- 

+ 1/3 + TCi 



(3.22) 



For arbitrary (3 and fixed scattering angle, we must solve for the relevant di- 
agonal basis of color structure, and determine the eigenvalues. At "absolute" 
threshold, /? = 0, where the heavy quarks are produced relatively at rest, we 
find 



r 
r 



5,11 
S,21 
5", 12 

S,22 



C-F 

TT 







Pis 

2^ 

TT 



J 



2^, 



2 



1 — ni — 



2^, 



(3.23) 



It is interesting to note that T^^^ is diagonalized |^ in a basis of singlet 
and octet exchange in the s channel. 



Csinglet 



Cl, 



1 ^C2 

Ci H . 

2N 2 ■ 



(3.24) 



with eigenvalues, Afj^J^^ = r§;7° and A^^^tet = ^3^2 ■ 

The general result in this s channel singlet-octet basis (including the 
contribution of the self-energy graphs, subtracting the z/*^*-* contribution as 
noted above and in a general axial gauge) is: 



.(1,8) 
5,11 

(1,8) 
5,21 

(1,8) 



-Cf[Lp + 1 + TTi + ln(2y^^)], 

2a. 



TT 



.ln(^ 



It, 



Ml 
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.(1,8) 
5,22 



IT 



4 In 



+ 



Ca 



-3 In 



-In 



tlUi 



+ L3 + 7Ti 



(3.25) 



Fs is also diagonalized in this singlet-octet basis for arbitrary j3 when the 
parton-parton cm. scattering angle is 6' = 90° (where ui = ti) with eigen- 
values that may be read off from Eq. ( p.25| ). 

It is of interest, of course, to compare the one- loop expansion of our 
results to known one-loop calculations, at the level of NLO. We give our 
result as a function of z, since the inverse transforms are trivial. The result 
is proportional to the Born cross section, the one-loop contributions from 
Qi^^^ and g2'^\ and r^^22'- In the DIS scheme the result is 



ff^QQ 



1 - z,ui,ti,s) 



(1) 



a. 



CTBo 



IT 1 — Z 



21n(l - z) + - 



+81n(^ 



2 - 2Re + 2 In I — 



+Ca 



-3 In 



+ ReL 



/3 



In 



(3.26) 

where we keep only terms that give rise to InA^. The logarithm of s//i^ 
describes the evolution of the parton distributions. This result for the cross 
section at fixed values of the QQ invariant mass cannot be compared directly 
to the one-loop results of because the latter is by contrast a single- 
particle cross section, where the singular behavior is given in terms of the 
variable S4, 

Si = {P2 + kf ~m^^2p2-k, (3.27) 

where k = pa + Pb — Pi — P2 is the momentum carried away by the gluon. At 
partonic threshold, both S4 and (1 — -2) vanish, but angular integrals over the 
gluon momentum reflect the difference in phase space between fixed S4 and 
fixed 1 — z ^ 2{pi + P2) ■ kj s. 

Nevertheless, the cross sections are kinematically equivalent in the /? — >■ 
limit, where we may make a direct comparison. Near /? = 0, our cross section 
becomes 



ff^QQ 



1 - z,ui,ti,s) 



(1) 



0"Bo 



1 



21n(l - z) 



3 
2 
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+2 In 



4m 



Ca ■ (3.28) 



S4, and this expression becomes 
. Even for /3 > 0, the two cross 



Near s = Am^, we may identify 2m?{l — z) 
identical to the /? — >^ hmit of eq. (30) of 
sections remain remarkably close, differing only at first nonleading logarithm 
in the abelian (Cf^) term, due to the interplay of angular integrals with 
leading singularities for the differing treatments of phase space. 

As for the Drell-Yan cross section, our analysis applies not only to abso- 
lute threshold for the production of the heavy quarks (/3 = 0), but also to 
partonic threshold for the production of moving heavy quarks. When (3 nears 
unity, however, the anomalous dimensions themselves develop (collinear) sin- 
gularities associated with the fragmentation of the heavy quarks, which in 
principle may be factored into nonperturbative fragmentation functions. 

Finally, it is easy to check that the anomalous dimension matrix for heavy 
outgoing quarks Eq. ( ^.19 ) reduces in the limit m — >■ to the anomalous 
dimension matrix for hght outgoing quarks, q{pa) + ^(Pb) q{Pi) + q{P2) , 
which is [p4| : 



S,ll 



a. 



In 



VabVl2 \ 



4 



\n{i'ai'b^i^2) + 2 — 27rz 



J^qq^qq 



T^qq^qq 
J- S,12 



-pqq^qq 
^ S,22 



qq^qq 
5,21 5 



— r 

2tt \VaiVb2J ' 

VabVl2 \ 
ValVh2) 



as 
2^ 



In 



71 



In 



2TTi 



V 4 



j - -ln(z/aZ/,,z/iZ/2) + 2 



-^^s^2 ■ 



(3.29) 



4 Applications to gg QQ 

In this section we give the results for the anomalous dimension matrix when 
the incoming partons are gluons and the outgoing partons are heavy quarks 



g{.Pa) + 9{Pb) Qi.Pi) + Q{.P2) 



(4.1) 



29 



For the sake of completeness we also give results for the case when the out- 
going quarks are light 



9{Pa) + gipb) q{p\) + q{P2) ■ 



(4.2) 



In Fig. 5 we show the UV divergent one-loop contributions to S for gg — > QQ. 

Our analysis is similar to the one in the previous section. We use the 
same integrals for the calculation of the oj'^^^ The eikonal rules for incoming 
gluons arc slightly modified and are given in Fig. 6. We choose the following 
basis for the color factors: 



Ci = 5 



ab r 



C2 = d 



abc 



rpc 
-'21) 



(4.3) 



Again, the counterterms for S are the ultraviolet divergent coefficients times 
our basis color tensors: 



51 — CiZs,ll + C2Zs,21 + C3Zs,31, 

52 — CiZs^u + C2^5,22 + ^3.^5,32, 
'S'3 = 13 + 02^5,23 + C3^S',33. 



(4.4) 



Our results for the contributions from vertex corrections in Fig. 5(a) to 
the anomalous dimension matrix when the outgoing quarks are heavy are: 



s,n — 



i 5,21 

r5,i2 

r5,22 



TT 



0, 

TT 

0, 

TT 



-Lr — Li 



+Ca 



In 



Vgb 

2 



L2 
1 



In 



Va2Vbl \ 
ValVb2) 



Cf(-L3 - Li - L2 + 1) + 



+ -\n{ValVb2Va2Vbl] 



Lp - In 



Ca 
2 

2m2 



In 



Vgb 

2 



- ln(z/aZ/fe) + 2-7Ti 
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(b) 



Fig. 5. One-loop corrections to Sji for heavy quark production through 
gluon fusion: (a) vertex corrections; (b) heavy quark self-energy graphs. 
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-iv7(-v-k+i8) -v7(v-k+i8) 




iv7(-vk+i8) v7(vk+i8) 

Fig. 6. Feynman rules for eikonal lines in gg QQ. The gluon momentum 
flows out of the eikonal lines. Group factors for vertices attached to incoming 
lines are those of three-gluon vertices, and for outgoing hues those of quark- 
gluon vertices. 



r5,23 



4 

r5,22 



S,31, 



(4.5) 



Again, the gauge dependence associated with the heavy quarks is eliminated 
once we include the self-energy graphs in Fig. 5(b). The contribution of the 
self-energy graphs (in the diagonal elements only) is, as before. 



a. 



TT 



Cf{Li + L2-2) 5ij. 



(4.6) 



In analogy to the previous section, we subtract i/'^'^/2, which gives an addi- 
tional —{as/7r)CA shift in the diagonal elements. In terms of the Mandelstam 
invariants s, ti, and Ui, the anomalous dimension matrix becomes (in a gen- 
eral axial gauge): 



r5,2i 



TT 

0, 



-CF{Lfs + l)-C,_ 



- In (4i/„i/fe) + m 
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r5,3i = 


as (ul\ 


rs,i2 = 


0, 


r5,22 = 


TT L 


r5,32 = 




4 AT 


r5,i3 = 




r5,23 = 


Ca 

5,31, 


r5,33 = 


r5,22- 



1) 



In (Auaiyb) - In ( 



(4.7) 



At threshold for the production of heavy quarks at rest, the anomalous 
dimension matrix becomes diagonal. In A'' = gauge (where z/^ = z/;, = 1/2) 
its eigenvalues are 



r 



5,11 
13^0 



S,22 



5,33 



TT 

as 

IT 



2/3 



Til 



(4.8) 



We also note that the matrix is diagonalized at 6' = 90° for arbitrary /5, with 
eigenvalues 



\i{e = 90°) 
X2{e = 90°) 

A3(^ = 90°) 



r5,ii( 

r5,22( 



90°) 
90°) 



+ 



cv 

^[~CF{Lf, + l)-CAT^l] 
TT 

71 

2 \ 



In 



t 



+ Lp- Tci) 



(4.9) 



Again, it is interesting to compare the one-loop expansion of our results to 
the one- loop calculations in In this case, our calculation is complicated 
by the fact that the color decomposition is not trivial as it was for qq. We 
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have to decompose the Born cross section into three terms according to our 
color tensor basis. After some algebra our result becomes 



X 



m^CpCA ln(l - z) + - In 



-2N^CfCa In 



+AN^Cf Cf - 



Ca 



ReLp + AN^Cl 



+8N'C'p ln(l - + - In — 



+CfCa{N^ -2) In 



2 



+2Cf Cf - 



where 



-2C|(iV2-l)} 



ti Ml Am?s 

Ml ti tiUi 



1 - 



rn?s 

tiUi 



(4.10) 



(4.11) 



and Kgg = {N"^ — 1)~^ is a color average factor. The logarithms of s/fi^ 
describe the evolution of the parton distributions. 

As we discussed in the previous section, our result cannot be compared 
directly to the one-loop results of p6| , but as /? — our expression becomes 



identical to the (3^0 limit of the sum of eqs. (36-38) in p6 



(1 - Z,Ui,ti, s)^^^|/3 



a: 



X 



Y^Kgg {m^CF{^CF-CA) 



ln(l - z) + - In 



Arrr 



+ {N^ + 2)CfCa~4:N^C],} . (4.12) 
Again, even for (3 > 0, the two cross sections remain remarkably close. 
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Finally, the anomalous dimension matrix for the case when the outgoing 
quarks are light is given by 



s,n 



5,21 

T^gg^Q^ 

^ 5,31 
5,12 

rgg^i^ 

5,22 



^ 5,32 



^ 5,13 



^ 5,23 
^ 5,33 



a, 
vr 

+Ca 
0, 

— In 



Cf 
In 



In 



VV2 

2 



1 



ln(z/iZ/2) + 1 



Til 



Vah 

2 



Til 



7T 

0, 

TT 



Va2Vbl \ 



ValVb2 



+Ca 



■1 



In 



Vl2 

2 



- - ln(z/iZ/2) + 1 



^ 5,31 ' 



Vab 
V12 



- - ln(z/aZ/fe) - ln2 + 1 



2 



-rgg-'m 
J- 5,31 5 



rgg^ii 

5,31 ' 



^ 5,22 ■ 



(4.13) 



We note that this matrix is also diagonalized at threshold and at ^ = 90°, 
and that the anomalous dimension matrix for heavy outgoing quarks, Eq. 
([4 .51) reduces to the anomalous dimension matrix for light outgoing quarks, 
Eq. ([4.13|) in the limit m — 0. 



5 Conclusions 

We have illustrated the application of a general method for resumming next- 
to-leading logarithms at the edge of phase space in QCD hard-scattering 
cross sections. We have given explicit results for singular distributions at 
partonic threshold in heavy quark production through light quark annihi- 
lation and gluon fusion, and for light quark production through gluon fu- 
sion. Possible extensions include, of course, dijet and multijet production. 
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and dijet or heavy quark transverse momentum distributions. Recent es- 
timates of the phenomenological importance of these nonleading terms to 
heavy quark production indicate that their contribution can be substantial 
[|1^ . In addition, an understanding of the all-orders perturbative cross section 
in QCD processes may shed light on power corrections in hadronic scattering 
[0, |2^, We hope that our method will be of use in improving the 

reliability of perturbative QCD calculations for hard scattering. 
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